Spatiotemporal chaos: the microscopic perspective 
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Extended nonequilibrium systems can be studied in the framework of field theory or from dy- 
namical systems perspective. Here we report numerical evidence that the sum of a well-defined 
number of instantaneous Lyapunov exponents for the complex Ginzburg-Landau equation is given 
by a simple function of the space average of the square of the macroscopic field. This relationship 
follows from an explicit formula for the time-dependent values of almost all the exponents. 



The search for connections between theories and quan- 
tities defined at different scales is the essence of statistical 
physics and the backbone of condensed matter theory. 
It nourishes the field of turbulence and spatiotemporal 
chaos where there is interest in finding connections be- 
tween dynamical characteristics such as fractal dimen- 
sions and Lyapunov exponents, and statistics of macro- 
scopic quantities such as correlation lengths. Such con- 
nections have not only a theoretical value but also impor- 
tant practical consequences because it is much easier to 
study macroscopic quantities than to obtain dynamical 
characteristics, especially in experiments [3 12]- 

In the last decade, statistical mechanics community 
has also been interested in relating dynamical charac- 
teristics of the system, e.g. Lyapunov exponents, KS en- 
tropy and fractal dimensions, with the macroscopic prop- 
erties, such as transport coefficients or entropy produc- 
tion, both in the classical and quantum systems [31-11 1 j . 
All deterministic systems studied within this perspective 
were finite dimensional. A natural question then arises if 
similar results can be obtained for spatially extended sys- 
tems. For instance, one would like to know the statistical 
properties of the fluctuations of phase space contraction 
rate and of the entropy production in driven fluid sys- 
tems. Infinite dimensionality of Navier-Stokes equations 
makes such inquiries a challenge, although some interest- 
ing conjectures have been proposed [T2"] . 

These considerations prompted us to consider the 
complex Ginzburg-Landau equation (CGL) and study 
the fluctuation properties of its phase space contraction 
and the connections to macroscopic quantities. CGL 
is a paradigmatic model of spatiotemporal chaos which 
in certain sense is intermediate between thermostated 
molecular dynamics models and realistic fluid systems. 
Due to its strong dissipative properties infinite dimen- 
sional CGL has a finite-dimensional attractor which can 
be appropriately described in terms of low spatial fre- 
quency Fourier modes [T5] . 

In this Letter we show that, even though the phase 
space contraction rate in CGL is infinite, one can con- 
sider contraction rate of volumes restricted to the inertial 



manifold, which is finite dimensional. This rate is equal 
to the sum of a finite number of instantaneous Lyapunov 
exponents. It turns out to be proportional to the macro- 
scopic mass of the field. Thus we have found out a direct 
relation between the "microscopic" sum of a finite num- 
ber of instantaneous Lyapunov exponents and "macro- 
scopic" mass of the field. We explore the structure of 
the spectrum of Lyapunov exponents and instantaneous 
Lyapunov exponents and show an approximate formula 
for large part of the spectrum of instantaneous Lyapunov 
exponents. The statistical properties of the fluctuations 
of phase space contraction rates and its relations to other 
macroscopic entropy-like quantities will be reported in a 
follow-up article |14j . 

We consider one-dimensional cubic complex Ginzburg- 
Landau equation on an interval of length L with periodic 
boundary conditions: 

A t = eA + (1 + t Cl )AA - (1 + ic 2 )\A\ 2 A , (1) 

where all the coefficients e,c±,C2 are real numbers. For 
convenience, let us restrict to a finite dimensional trun- 
cation in Fourier base with N = 2K modes and write 

K 

A(x,t) = A n (t)e t2mix/L . 

n=-K 

From eq.Q we obtain 

: / 27m \ 2 , 

A n = £A n-[-FT) (l + ici)A n 

-(l + ica) J2 A ^A* m . (2) 

fe-f-Z— m=n 

Note that Ak = A-k due to periodicity. Writing A n — 
B n +iC n where B n and C n are real we derive a formula for 
the phase space contraction rate a — div^ A — J2 n Wf~ + 

as well as the normalized phase space contraction 
rate a := a/N mo ^ es , where N mo ^ es = 2N = AK is the 
number of real modes under examination: 

a , v [2tt\ 2 N 2 + 1 
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where (g) = 
L/N we get 



(l/L) dx\A\ 2 — J2k \Ak\ 2 - Using a = ume spanned by Ui- Its contraction rate is 



a = e-2(g)-^(l 
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The beauty of this result connecting the average macro- 
scopic field (g) to the microscopic normalized phase space 
contraction rate a is jeopardized by the last term that di- 
verges when the spatial resolution N is increases. How- 
ever, increasing the resolution only adds high frequency 
modes which are strongly damped. Wc show below that 
their contribution can be isolated and removed, as it is 
the case for zero-temperature entropy in spin systems. 

We conjecture that that there is a distinguished dimen- 
sion such that the contraction rate of volumes restricted 
to this dimension are always finite and connected to the 
space averaged g in a simple manner. These volumes are 
defined by the sum of an appropriate number of instan- 
taneous Lyapunov exponents. Before supporting these 
claims let us recall the definitions of Lyapunov exponents 
and instantaneous Lyapunov exponents, and show how 
they connect to the volume contraction rates. 

Consider a continuous time dynamical system defined 
by a set of differential equations x = F(x),x € M. n . The 
solution of the system is given by the flow xt = $*(xo), 
t £ 1. Then the growth of an infinitesimal perturba- 
tion 5xq around xq is governed by the linearization of 
the flow Sx t = Aco^* ' $ x o = M(t,xo) ■ 5xq- The fun- 
damental matrix M(t,xo) governing this growth is the 
solution of the equation M(t,xo) = J(t,Xo) • M(t, xq), 
where J(t,xo) — ^(^'(^o)) is the Jacobi matrix of 
partial derivatives of the field velocity. Oseledec ma- 
trix (M(xo, £)^M{xq, t)) 1 / 21 has n positive eigenvalues 
Ai(xo,t) which we order by size A x > A 2 > . . . > A n . 
Lyapunov exponents Aj(xo) are defined as logarithms 
of eigenvalues of long time limit of Oseledec matrix 
Ai(xo) := limt-j.oo In A(xo, t). For an ergodic system, Lya- 
punov exponents are the same for almost every initial 

point nunsi. 

To define instantaneous Lyapunov exponents [17] Hi 
consider volume Vfc(t) of a parallelogram ui(xo,t) A 
U2(xq, t) A ... A Ufe(xo, i), spanned initially by k orthogo- 
nal vectors Hi attached at xo, travelling along the trajec- 
tory; Ui € R n . Its evolution is given by the fundamental 
matrix, i.e. «j(xo,t) = M(xo,t)Ui. Then the k x n ma- 
trix U = [tii, ... , Uk] can be uniquely decomposed into a 
product of k x n orthogonal matrix Q and upper-diagonal 
k x k matrix R (QR decomposition) 



U = QR= [Qi, 
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We define instantaneous exponents by jJL k {t) := Cfc(i) — 
Cfc-i(i). They depend on the initial point and on the 
initial vectors Ui. However, for almost all initial vectors, 
the first vector with time aligns along the most unstable 
direction, the first two vectors span the fastest stretching 
2d volumes, and so on. Therefore, after some time the 
vectors become almost independ of the initial directions 
modulo degeneracy, and consequently the instantaneous 
Lyapunov exponents characterize the trajectory. 

In practice, we propagate the vectors by finite time 
steps at each time reorthogonalizing the set. Thus start- 
ing from Qo = U we move to U\ = M(dt)Q Q = Q\R\. 
Then U n+ i = M(dt,x(ndt,x ))Q n = Qn+\R n +i- Thus 
we have R(n-dt) = R n -. . .-Ri and fik{n-dt) sw 4r \n[R n ]kk- 
Time averages of fa are sorted in decreasing order and 
equal to the usual Lyapunov exponents Xi [HI [16] . 

To estimate the values of the instantaneous Lyapunov 
exponents, we consider an initial perturbation SA n tan- 
gent to a single mode . Inserting A n = A^ + SA n into 
eq. 0, we obtain 

dSAJdt = (e - (1 + ici)q 2 - 2(1 + ic 2 ){g))dA n (4) 
+(1 + ic 2 )(a n + iP n )A* n + f(SA))} (5) 

where q := 2irn/L, f(6A) is a linear function of {6 A} not 
depending on 5A n or <L4* , and a n , j3 n € R stand for the 
real and imaginary part of the time dependent sum 
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Rewriting the equation for real and imaginary parts of 
A n = B n + iC n , we can obtain short time evolution of 
tangent vectors 



SB n (dt) 
SC n (dt) 



[a Q L + a. t <7i] 



5B n (0) 
SC n (0) 



where <7j are the Pauli matrices |18j and oq — 1 + (e — 
q 2 -2{g))dt, a x = P n +c 2 a ni a y = i(c x q 2 + 2c 2 (g)), a z = 
ot n — C2j3 n . Then the eigenvalues of M'M(dt) are A± = 
1 + 2(e -q 2 - 2(g)) ± yTT^|| Ef=J(k|„|) A n-^ n+J \ 
which gives extremum possible values of instantaneous 
Lyapunov exponents 



/i„± =e-q 2 -2(g)± Jl + 



K-\n\ 



An-jA n 



j=-(K-\n\) 



The product of the diagonal elements of R gives the vol- 



Observed values depend on initial vector [SB n SC n ] T and 
are between /i n ±. To find out what is the contraction 
rate of the 2d volumes in SA n plane consider the action of 
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M(dt) restricted to 5A n on a pair of initially orthogonal 
vectors [vi, V2] in this plane. The volume of M(dt)[v\,V2\ 
is given by the determinant, and since det[ui, 1*2] = 1, we 
have 

\{pni+Hva) = Um —Indet M(di) = e-q 2 -2(g). (6) 
z dt— >o at 

Therefore, at any time we predict that the sum of the two 
instantaneous Lyapunov exponents for perturbations in 
the plane tangent to any Fourier mode should be given 
by above formula. 

It is not a priori obvious that this prediction holds for 
any exponents calculated for volumes evolved over long 
time span, since we have considered the evolution along 
(A n ,An) only. In fact, realistic evolution mixes all the 
modes via the nonlinear term in ([I]) , and vectors align ar- 
bitrarily in phase space along their evolution. However, 
our numerical simulations show that there is only a finite 
number of modes W, which we call "active", behaving 
in apparently random (though smooth) fashion, which 
disobey the above prediction (Figure The remain- 
ing exponents in their time course oscillate around Q 
and at every instant the sum of the four modes for ±n 
is given by Sorting the instantaneous Lyapunov ex- 
ponents according to their time averages (i.e. Lyapunov 
exponents Aj), the "active" instantaneous Lyapunov ex- 
ponents are the first W curves. The "active" modes in- 
clude those tangent to the inertial manifold, the remain- 
ing exponents describe the decay towards the attractor. 

There remains the question of obtaining W. Our 
numerical results show that W is the smallest number 
of the form An + 2 greater or equal to L, i.e. W = 
2 + 4|~(L — 2)/4], where [~n] is n rounded upwards. This 
is larger than the Kaplan- Yorke dimension of the attrac- 
tor, and sometimes much larger. Therefore we believe 
our procedure probe the fluctuation of volumes on the 
inertial manifold, not on the attractor, which is a subset 
of it. W correspond in real space to a given size, coher- 
ent objects of size larger then L(W — 2)/(4.2n) drives 
the dynamics while structures of lower size are slaved to 
them. 

We have studied the spectra of Lyapunov exponents 
and instantaneous Lyapunov exponents for CGLe for sev- 
eral parameter values and truncations to 32, 64, 128 
and 256 Fourier modes or equivalent numbers of spatial 
points. Figure [T] shows the time dependence of the first 
46 instantaneous Lyapunov exponents computed along a 
trajectory for CGLe with L = 10n (q = 0.2), cj = 4, 
C2 = —4 [TJ5] for a computation with 128 Fourier modes; 
both groups of exponents are clearly visible. There is 
a constant difference between theoretical prediction and 
the numerical value of the order 0.15 for the first inactive 
exponents (n = 9 in this case), dropping to around 0.01 
for n = 25. 

Figure [2] shows the spectrum of average instantaneous 
Lyapunov exponents, i. e. of Lyapunov exponents, for 
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FIG. 1: Time dependence of all the 34 "active" and the largest 
12 "inactive" instantaneous Lyapunov exponents for CGLe 
with L — Wn, c\ — 4, C2 = —4. The bottom figures show two 
subsets of curves from the upper plot, the first 34 (left) and 
the next 12 exponents (right). Prediction (JsJ) is also plotted. 

CGLe with the same parameter values. Crosses are theo- 
retically predicted values ([6]), circles are numerical values. 
The staircase structure in the spectrum is well approxi- 
mated by (|6j). 

To separate the changes of volumes on the inertial 
manifold from the trivial contraction of infinite dimen- 
sional phase space onto finite dimensional inertial man- 
ifold we consider contraction of W dimensional volumes 
given by the sum of all the nontrivial instantaneous Lya- 
punov exponents. Since the sum of all the instantaneous 
Lyapunov exponents in the Galerkin representation is 
equal to the phase space contraction rate ([3]), and since 
the "inactive" exponents on the average follow the aver- 
age field Q, the relevant value is 

1 , , 7T 2 (W 2 -A) 

Active := = 2^ Mi = e - 2 (e> ^ ■ ( 7 ) 

This is approximately <7 ac ti V c ~ £ — 2(g) 12L2 ~ 

—2(g). Figure [3] compares the evolution of 0.15 + 

W 2^=1 Mi with e - 2(g) \ 2L2 ' . 

To summarize, we have shown that the phase space 
contraction rate in Galerkin approximation to the com- 
plex Ginzburg-Landau equation is given by a simple func- 
tion of the spatial average of the squared modulus g of 
the solution. The divergence occuring when increasing 
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FIG. 2: Spectrum of Lyapunov exponents for CGLe with L = 
IOtv, ci = 4, C2 = —4. iV = 64 modes, there are 128 exponents 
(all shown in the inset). The quadruplet structure of the 
lower exponents comes from the real and complex parts of 
corresponding positive and negative Fourier modes. Crosses 
are theoretical prediction circles are numerical values. 
There are W = 34 "active" exponents. 
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FIG. 3: Numerical test of formula ffb. The time evolution of 
the sum of the first W=34 instantaneous Lyapunov exponents 
(active modes; in red) are compared with the time evolution 

of e — 2(g) — n ^ L 2~ 4 " > (in blue). The second term has been 
shifted by -0.15. 



spatial resolution can be removed by restricting the con- 
tracting volumes to a finite number of dimensions. The 
corresponding volume contraction rate is given by the 
sum of a finite number of Lyapunov exponents which 
time-behavior is non-trivial. We have identified a natu- 
ral division of the spectrum into the part corresponding 
to the dynamics on the inertial manifold and the other 
part corresponding to the modes decaying towards the 



attractor. Instantaneous Lyapunov exponents in the sec- 
ond part are approximately given by simple functions of 
the square of the Ginzburg-Landau field Q. We have 
found a formula for the volume space contraction rate 
on the inertial manifold The formula bridges the 

gap between the dynamical systems picture of the CGLE 
(volumes contracting in the phase space and instanta- 
neous Lyapunov exponents) and the macroscopic picture 
(spatio-temporal solution) . 
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